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Abstract 

In this paper, the problem of rate control in multi-user multiple-input multiple-output (MIMO) 
interference systems is formulated as a multicriteria optimization (MCO) problem. The Pareto rate 
region of the MCO problem is first characterized. It is shown that for the convexity of the Pareto 
rate region it is sufficient that the interference-plus-noise covariance matrices of multiple users 
with conflicting objectives approach identity matrix. The latter can be achieved by using either 
orthogonal signaHng, tame-sharing, or interference cancellation strategies. It is argued that, in the 
case of high interference, the interference cancellation is preferable in order to increase the Pareto 
boundary and guarantee the convexity of the Pareto rate region. Secondly, the Nash bargaining (NB) 
is applied to transform the MCO problem into a single-objective problem. The characteristics of 
the NB over MIMO interference systems such as the uniqueness, the existence of the NB solution 
and the feasibility of the NB set are investigated. It is found that when the NB solution exists, 
the sufficient condition for the corresponding single-objective problem to have a unique solution is 
also that the interference-plus-noise covariance matrices of multiple users approach identity matrix. 
A simple multi-stage interference cancellation scheme, which leads to a larger convex Pareto rate 
region and, correspondingly, a unique NB solution with larger user rates compared to the orthogonal 
and time-sharing signaling schemes, is proposed. The convexity of the rate region, the effectiveness 
of the proposed interference cancellation technique, and the existence of the NB solution for MIMO 
interference systems are examined by means of numerical studies. The fairness of the NB solution 
is also demonstrated. Finally, the special cases of multi-input single-output (MISO) and single-input 
single-output (SISO) interference systems are also considered. 

Index Terms 

Multicriteria optimization, Pareto region, Nash biffgaining, game theory, MIMO interference 

channel. 
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I. Introduction 

Multicriteria optimization (MCO), which is featured by the need to simultaneously op- 
timize multiple noncommensurable or even incompatible objectives, has been widely used 
in automatic control, telecommunications, economics, and many other fields of engineering 
and science [1]. An MCO problem usually admits infinite number of noninferior solutions, 
which form a Pareto boundary and encompass a Pareto region [2]. A noninferior solution in 
the Pareto boundary is considered to be Pareto optimal in the sense that no other solution 
can improve the performance of some objectives without deteriorating other objective(s). The 
characteristics of the Pareto region are of particular interest in terms of shedding light on 
solving the corresponding MCO problem. 

Two main processes are involved for solving the MCO problem, namely, computing and 
decision making [3]. In the former process, a set of noninferior solutions for the MCO 
problem is obtained, while, in the latter one, a final solution among the noninferior solutions 
is determined according to certain preferences of the system under optimization. Depending 
on the principle according to how these two processes interact with each other in search for 
a compromise solution, the MCO problem can be solved using one of the following three 
approaches [3]: priori, posteriori and progressive articulation of preferences. As the name 
suggests, the decision maker expresses its preferences prior to the computing process in the 
first approach. These preferences are often given in the form of an utility function aggregating 
all objectives. As for the posteriori case, the preferences of the decision maker are not known 
before computing. Thus, all the noninferior solutions for the MCO problem should be found 
before decision making. Unlike the first two approaches, in the progressive articulation of 
preferences, computing and decision making are performed in an alternating manner. In each 
iteration, a population of noninferior solutions is maintained during computing; then, the 
decision maker chooses several suitable compromise solutions for the computing process 
in next iteration to produce a new generation of noninferior solutions. After rounds of 
iterations, the final solution is obtained. Evolutionary algorithms are usually employed in 
this approach [4]. 

Among the three abovementioned approaches for solving the MCO problem, the priori 
articulation of preferences is most commonly used due to its simplicity and effectiveness. It 

can be further broken down into three broad classes of methods: scalarization of multiple 
objectives - the MCO problem is transformed into a single-objective optimization problem 
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by scalarizing the multiple objectives into a single- valued utility function [3]; prioritization 
of multiple objectives - a certain priority is assigned to each objective and the problem is 
optimized according to the priority order of its objectives [5]; goal-value method - a goal value 
(usually the minimum attained performance) is assigned to each objective before optimization 
to indicate its corresponding desired level of performance [6]. 

The rate control problem in multi-user interference systems can be formulated as an 
MCO problem. The convexity of the Pareto rate region of the MCO problem is a desirable 
feature for the interference system. Various approaches have been applied to convexify a 
rate regiorQ [7] [8]. These approaches can be broadly divided into two categories. The first 
one is based on the use of a convex hull (see, for example [7]), which corresponds to 
the so-called time-sharing strategies. However, the time-sharing rate boundary is achievable 
only in terms of average rate rather than instantaneous rate. The other approach is based 
on using orthogonal stategies (signaling schemes), e.g., time division multiplexing (TDM) 
and frequency division multiplexing (FDM) [8]. Orthogonal signaling schemes significantly 
simplify the MCO problem, but they usually lead to smaller achievable rate regions. The 
corresponding rate loss can be significant when the interference is low. Therefore, it is 
desirable to analyze the convexity of the true Pareto rate region only employing pure non- 
orthogonal strategies. 

To solve the MCO problem in multi-user interference systems, several methods have 
been developed in the literature based on the priori articulation of preferences. The most 
representative one is sum-rate maximization, e.g., the sum-rate maximization for multi-user 
multiple-input multiple-output (MIMO) interference systems in [9] . The obvious drawback of 
this method is the lack of fairness, since the performance of users with bad channel conditions 
are always sacrificed. The fairness can be improved by weighted sum rate approaches like 
in [10], but it is not always easy to determine the weighted coefficients for all users. The 
preference could also be based on proportional fairness like in [11], where the fairness of a 
user is proportional to its channel conditions. Recently, game theory has been increasingly 
used in the control and optimization for communication systems like in [12]. Specifically, 
the Nash bargaining (NB) from cooperative game theory has been widely used to solve the 
MCO problem for multi-user interference systems. NB is an effective scheme to balance the 



'when referring to rate region/boundary, we mean Pareto rate region/boundary hereafter, if there is no particular other 
clarification. 
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fairness of individual users and the system-level performance [13]. Representative examples 
in the literature include [8] and [14], for single-input single-output (SISO), [7] and [15] for 
multiple-input single-output (MISO), and [16] for MIMO interference systems. In [16], a 
practical suboptimal algorithm for finding the NB solution was designed by exploring the 
gradient projection method [9]. However, little research has been done to characterize the 
pure strategy based NB for MIMO interference systems. 

In this paper, we apply MCO to rate control problem in multi-user interference systems. 
Specifically, we formulate the cooperative rate control of MIMO interference systems into 
an MCO problem. Firstly, the Pareto rate region of the MCO problem is characterized. It is 
proved that the interference-plus-noise covariance matrices approaching the identity matrix is 
a sufficient condition for the convexity of the rate region. Moreover, a significant implication 
is found that when interference is high, interference cancellation techniques are preferable 
in making the rate region convex. Secondly, the MCO problem for rate control in multi-user 
MIMO interference systems is converted to a single- objective Nash-product maximization 
problem using NB. The characteristics of the NB over MIMO interference systems such as 
the uniqueness of the NB solution, the feasibility of the NB set and the existence of the NB 
solution are studied. Furthermore, we propose a multi-stage interference cancellation scheme, 
which leads to a larger rate region and, correspondingly, an NB solution with larger rates 
compared to orthogonal and time-sharing signaling based approaches. More importantly, 
it also guarantees the convexity of the rate region when the interference is high, which 
consequently admits a unique NB solution. 

The remainder of this paper is organized as follows. The MCO problem and Nash-product 
maximization problem for multi-user MIMO interference systems are formulated in Section II. 
The main results including the characterization of both the Pareto rate region and the NB 
solution, and the elaboration for the proposed multi-user interference cancellation technique 
are given in Section III. The convexity of the rate region, effectiveness of the proposed 
interference cancellation technique, fairness, existence and uniqueness of the NB solution for 
MIMO interference systems are exemplified in Section IV via numerical studies. In Section 
V, some special cases are reviewed. Conclusions are drawn in Section VI. 
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II. Problem Formulation 
A. MCO in MIMO Interference Systems 

Consider an M-user MIMO interference system in which all users transmit on the same 
wireless channel simultaneously. The transmitter and receiver for each user are equipped with 
Nt and A^^ antennas, respectively. The iV^ x 1 complex baseband signal vector received by 
user i (i = 1, 2, . . . , M) can be written as [17] 

M 

Yi = y/plHiiXi + ^ y/r]~iiijyij + (1) 

where pi is the normalized signal-to-noise ratio (SNR) for user i; rjij — 1, - • • , M, i j) 
is the normalized interference-to-noise ratio (INR) from transmitter j to receiver i; Ha and 
H.y are the A^^ x Nt channel matrices from transmitter i and transmitter j to receiver i, 
respectively; Xj is the A^^ x 1 transmitted signal vector for user i; and is the iV^ x 1 
independent and identically distributed (i.i.d.) additive white Gaussian noise (AWGN) vector 
of user i with zero mean and covariance matrix E[ninf] — I. Here, E[-\ stands for the 
expectation operator, (•)^ denotes the Hermitian transpose operation, and I is an A^,. x A^,. 
identity matrix. 

We assume that: (i) each transmitter/receiver transmits/receives symbols independently; (ii) 
the co-channel interference from other users is unknown and treated as noise, i.e., no interfer- 
ence cancellation techniques are employed by receivers; (iii) the channel varies sufficiently 
slowly and is considered as time invariant during the period of each symbol transmission. 

The mutual information for user i can be expressed as [18] 

/,(Q) = log2det(I + p,HnQ^H,^R=•), ^ = 1,...,M (2) 

where Qj = £'[xjxf ] is the Hermitian positive semi-definite (PSD) transmit covariance matrix 
of the input signal vector for user i, i.e., Qj >z 0, and 

M 

R_, = I+ %H,,Q,Hg, i = l,...,M (3) 

is the interference-plus-noise covariance matrix for user i. We define Q = (Qi,...,Ql) as 
a set of transmit covariance matrices. Since the transmission of each user is power limited, 
the following trace constraint applies to Qj 

tr(Q,)<Pi, i = l,...,M (4) 
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where tr(-) denotes the trace operator. We also assume that each transmitter i has the full 
knowledge of the channel, SNRs, INRs and transmit covariance matrices of all the other 
transmitters. 

The rate control objective in the MIMO interference system is to maximize the rateElof 
each user by optimizing their transmit covariance matrices Qj {i = 1, . . . , M) under the 
trace constraints given by dH). Therefore, rate control in the multi-user MIMO interference 
system can be formulated as the following MCO problem 

max/,(Q) i = l,...,M 

Q 

subject to Qj ^ 0, i = 1, . . . , M 

tT{Q^)<p^, i = l,...,M. (5) 

B. Scalarization of the MCO using NB 

According to game theory, a game consists of three elements: players, strategy and util- 
ity [19]. Players are rational parties involved in the game. Strategy stands for actions or 
behaviors taken by players. Utility is usually defined in the form of a certain performance 
metric for players. The MIMO interference system delineated above can be modeled as a 
MIMO interference game, whose players are the users in the MIMO system. The rate of 
each user represents the utility of the corresponding player. The transmit covariance matrix 
of each user forms the strategy space of each player. 

A game can be classified as either competitive or cooperative according to the cooperation 
scheme among players. In a competitive game, as the name suggests, all the players compete 
with each other selfishly and rationally. Players neither communicate nor cooperate with each 
other during the game. A steady state in a competitive game for which each player can not 
improve its utility by unilaterally changing its own strategy is called the NE [19]. For a 
MIMO interference game, the NE can be mathematically expressed as 

= 1, . . . , M, Q, ^ 0, tr(Q,) < : J,(Q*, Q%) > /,(Q„ Q*_J (6) 

where Q* and Q^i^ denote the transmit covariance matrices of the NE for user i and for all 
the other users except i, respectively. The transmit covariance matrix of each player leading 
to the NE can be found via iterative water filling (IWF) as [20] 

Q* = U,(/iJ-Dri)+Uf, z = l,...,M (7) 



^Hereafter, when referring to rate, we mean mutual information obtained using 
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where UjDjUf = piH^^HzlHu is the eigendecomposition of piH^^HzlUa, Ui {i = 
1, . . . , M) is the unitary matrix of eigenvectors, Dj (i = 1, . . . , M) is a diagonal matrix 
of eigenvalues, and fii denotes the power level given by IWF. 

Generally, the NE is not optimal from the system point of view due to its competitive 
and selfish nature. Whereas, in a cooperative game, all the players negotiate with each other 
prior to the game, which usually results in utility improvement [21] [22]. There exist many 
cooperative game-theoretic approaches. In this paper, we restrict our attention to the NB as 
the simplest and also the most popular one [13] to scalarize the multiple objectives of the 
MCO problem ©. 

The NB is well defined in a convex rate region [13]. In the context of the MIMO 
interference game, for the case when the rate region is convex, the bargaining set, i.e., the 
set of available strategies for user i, can be expressed as 

5' = {Q.|Q. ^0, tr(Q,)<p„ and/,(Q)>/f^,2 = l,...,M} (8) 

where Lj^^ is the utility of the NE for user i. It also has a disagreement point, which is 
defined as the state that players resort to when the cooperation fails. Usually, the NE is taken 
as the disagreement point in NB. By applying the NB, the multiple objectives in the MCO 
problem © can be scalarized. Then, the MCO problem can be transformed into the following 
single-objective optimization problem 

M 

maxn(/.(Q)-/r) 

1=1 

subject to Qi ^ 0, i = 1, . . . , M 

tr(Q.)<p„ z = l,...,M 

/i(Q)>/f^, t = l,...,M. (9) 

It is worth noting that the NB corresponds to proportional fairness and the Nash product 
]^*f^(/j(Q) — /j^^) is converted to the rate product of all users when Ij^^ = for all 
i. An intuitive explanation of the above optimization problem is that the NB introduces 
a cooperation scheme among all MIMO users and regulates their transmissions (transmit 
covariance matrices) under the power constraint. On one hand, it guarantees that the utility 
of each MIMO user is not less than the one given by the NE. On the other hand, it maximizes 
the Nash product of the whole MIMO system. Therefore, it provides a good trade off between 
the fairness requirements to individual users and the overall performance of the whole system. 
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In the sequel, we first characterize the Pareto rate region of the MIMO interference systems. 
Then, we study several characteristics of the NB in MIMO interference systems such as 
the uniqueness of the NB solution, the feasibility of the NB set and the existence of the 
NB solution. Finally, a simple multi-stage interference cancellation scheme is proposed to 
convexify the rate region. 

III. Main Results 
A. Characterization of Pareto Rate Region 

The convexity of the rate region is a desirable feature for communication systems. Its 
immediate merit is that it usually yields larger rate region compared to a nonconvex one. 
More importantly, the convexity of the rate region is a necessary condition to ensure the 
convexity of the MCO problem. It simplifies the MCO problem and makes it solvable and 
practically realizable. There exist two broad approaches to convexify a rate region. The most 
common and direct approach is based on the use of a convex hull, which corresponds to the 
time-sharing signaling. In this case, all users agree to use a set of transmitting strategies for 
certain fraction of time and to use another set of strategies during the rest of the time (see, 
for example [7]). However, time-sharing signaling leads to a convex rate region in terms of 
the average rate (over time). Thus, the time-sharing rate boundary may not be achievable in 
terms of instantaneous rate [7]. 

Another approach is based on the use of orthogonal signaling among cooperative users. 
The examples of orthogonal signaling are and FDM. In this case, users agree to split 

the degrees of freedom (time or bandwidth) into several orthogonal parts and each user uses 
one part only. The orthogonal signaling significantly simplifies the MCO problem, but the 
associated rate loss may be significant, especially when interference is low. Instead of using 
either time-sharing or orthogonal signaling, we investigate the convexity of the true Pareto 
rate region only using pure strategies. The true Pareto rate region is both achievable in terms 
of instantaneous rate and usually larger than the orthogonal signaling-based counterpart. In 
the context of the MIMO interference systems, we derive the following sufficient condition 
ensuring the convexity of its rate region. 

Proposition! : If the interference-plus-noise covariance matrices R_i — > I (z = 1, ■ ■ ■ , M), 
the Pareto rate region of the MIMO interference system is convex. 

^TDM distinguishes itself from the time-sharing signaling/convex hull in the sense that for the 2-user case the TDM rate 
boundary is a line connecting two single-user points, while time-sharing signaling corresponds to the outmost lines between 
points in the Pareto rate boundary as shown in Fig. 2. 
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Proof: It is straightforward to prove that the rate region is convex when the utility 
functions of all users /^(Q) (i = 1, ■ ■ ■ , M) in Q are concave functions of Q. Let us prove 
the concavity of the function /^(Q), i = 1, . . . , M under the condition that interference-plus - 
noise covariance matrices approach I. Note that a function g(x) is concave if and only if (i) 
/(t) = g{txi + (1 — t)x2), < t < 1 is a concave function of t for any feasible Xi and X2, 
which is equivalent to f"{t) = d?f{t)/dt^ < 0; and (ii) the domain of g{x) is convex [23]. 

We consider the following convex combination of two different sets of transmit covariance 
matrices: (Xi, . . . , X^) and (Zi, . . . , Z^), that is 

Q(t) = (l-t)(Xi,...,XM)+t(Zi,...,ZM) 

= (Xi, . . . , Xjv/) + t(Zi — Xi, . . . , 7am — Xm) 

= (Xi,...,XM) + t(Yi,...,YM) (10) 

where < t < 1 and Yj = Zj — Xj (i = 1, . . . , M). Let us expand the utility function Q as 

m) = mit)) 



log2det (l + p,H,,Q,HfR-i) 

1 ^^ det(R_, + p,HiiQiHg) 
ln2 ^ det(R_i) 



1 M. 



Applying the well known property of matrix differential calculus [24], i.e., 

4- Indet (A(x)) = tr f A{x)-'^^ 
ax \ ax 

we can obtain the first derivative of as 



(11) 



(12) 



1 



^ In det(R_i + piHiiQiHf ) - ^ In det(R_, 



tr (R_,+piH,,Q,Hf) 



(iR_j 
dt 



+PjHjjYjH, 



H 



-tr r: 



-iC^R-i 



.(13) 



In ([T2l) . A(x) is a matrix function of scalar parameter x. Using ([T3l) and applying two other 
properties of matrix differential calculus [24], i.e., 

fdA{x) 



-^tr (A(x)) = tr 
dx 



-A(x' 



\ dx 
_^dA{x) 



dx 



Aix)-\ 



(14) 
(15) 



the second derivative of fi{t) can be expressed as 



1 

h[2 



tr (-(R_. + M,)~^ (^ + ^') 



M, 



dR_ 
~dt 



+ tr(R:l^R-i^^ 



dt 



dt 



(16) 
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where M, = piHiiQ^Hf and N, = p^HuYiHfl . 

Let Aj = (R_j + Mj)~^ and Bj = dK^i/dt + Nj. Since Aj ^ 0, there exists a matrix Cj 
such that Aj = CjCf^ . Thus, the first trace in the right hand side of (fT6l) can be written as 

tr (-AjBjAjBi) = -tr (C,Cf B.QCf B,) 

= -tr (Cf B,C,Cf B,Q) 

= -tr ((Cf B,Q) (Cf B,;Q)'') < (17) 

The last equality in (flTI) is obtained using the fact that Bj is Hermitian. The inequality in 
^ holds due to the fact that (Cf BjCj) (Cf BjCj)^ >z 0. 

When R_j — )■ I, in the right hand side of (fT6l) . dH^i/dt = Vij'^ij^ij^fj approaches 
0. Then, in the first trace of (fT6l) . (R_j + Mj)^^ and {dH-i/dt + Nj) are dominated by 
(I + Mj)"^ and Nj, respectively. Note also that the second trace can be ignored as compared 
to the first one. Therefore, when R_j — )• I, f"{t) < 0. 

The domain of the utility function /i(Q) for the MIMO interference system is {Qj|Qj ^ 
0, tr(Qj) — Pi < 0, i = 1, . . . , M}, which is clearly convex. Therefore, the utility function 
/j(Q) is concave and consequently the rate region is convex when R_j — )• I. □ 

The condition that R_j — )• I means that Vij^ij^ij^ij ~^ 0- This corresponds to the 
condition that INRs r]ij (i, j = 1, ■ ■ ■ , M, i ^ j) in ([T]) are sufficiently small. Hence, the 
immediate implication of Proposition 1 is that the rate region is convex when the interference 
is low. 

In summary, we know that the orthogonal signaling is a sufficient condition for convexity of 
the rate region [8]. At the same time. Proposition 1 suggests that strict signal orthogonality 
is not a necessary condition for convexity of the rate region. Another sufficient condition 
which guarantees the convexity of the rate region is R_j — )■ I (i = 1, ■ ■ ■ , M). It generalizes 
the requirements of the time-sharing and orthogonal signaling, when constructing a convex 
rate region. More importantly, the following remark of a practical significance can be made. 

Remark 1: When the interference is high, interference cancellation techniques outperform 
orthogonal signaling techniques for convexifying a rate region in the sense that they lead to 
a convex rate region of a larger size. 

Orthogonal signaling is a simple and widely-used method to produce a convex rate region 
at the cost of rate loss. When the interference is high, applying interference cancellation 
techniques to a MIMO interference system eventually transforms a high-interference system 
into a low-interference system, which consequently leads to a convex rate region according 
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to the Proposition 1 and results in a larger rate region compared to the true Pareto rate region. 
Therefore, interference cancellation techniques outperform orthogonal signaling techniques 
in terms of the size of the resultant rate region. 

Various interference cancellation techniques [25] can be applicable to interference sys- 
tems. The widely used multi-user detection technique [27], including sucessive interference 
cancellation (SIC) and parallel interference cancellation (PIC), is one suitable candidate for 
interference systems. Specifically, for MIMO interference systems, beamforming is another 
effective interference rejection technique. Transmit and receive beamforming can be per- 
formed in MIMO transmitters and receivers to mitigate interference, respectively. However, 
providing the interference cancellation capabilities for the interference system requires in- 
crease of both computational complexity and communication overhead among users. Take 
SIC and PIC for instance, these two approaches need to decode and reconstruct the interfering 
signals before substracting them from the received signal. The decoding and reconstruction 
processes definitely complicate the intereference system. Besides, channel information for the 
interfering channels is also required by the decoding and reconstruction. Moreover, additional 
communication overhead among receivers is necessary as well when receiver cooperation is 
employed to further enhance the interference cancellation. The additional complexity and 
communication overhead are paid off by the increase of rate region. 

B. Characterization of NB Solution 
1) Uniqueness of NB Solution: 

It is straightforwad to see that the optimization problem given by ^ is identical to 
max In ||ln^ 2 J] (/,(Q) - Jf^) j 

subject to Qj ^ 0, i = 1, . . . , M 

tr(Q.)-pi<0, ^ = 1,...,M 

Jf^-/,(Q)<0, t = l,...,M. (18) 

Note that ([TSl) is a convex optimization problem if and only if its objective function is concave 
and its constraint set is convex [23]. 

Proposition 2 : When the interference-plus-noise covariance matrices R_j — t- I (i = 
1, ■ ■ ■ , M) , the optimization problem ([TSl) is a convex problem. 
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Proof: First, by using the same methodology as the one used to prove Proposition 1, let 
us show that the objective function in (fTSl) is concave under the condition that the interference- 
plus-noise covariance matrices R_j approach I, Vi Considering the convex combination in 
(flOl) . the objective function of (fTSl) can be expanded using (O as 

/ M 

/(t)=ln(ln^2n(/.(QW)-/r^ 

i=l 



In (hi 

1=1 ^ 



det(R_j + PiHiiQiHf ; orA^S 



det(R_i 



-ln2^f^ . (19) 



Let us define det(R_, + aH,,Q,H^) , 

Ti = In ^ ' , - In 2Jf^. (20) 

det(R_j) 

Using (fT2l) . (fT4l) and (fTSl) . we obtain the second derivative of f{t) as 

A/ 

/"(t) = 5^a + /3 + 7 (21) 



i=l 



where 



tr (- (R_, + Mi)-^ f + N,) (R_,; + M,)-^ + N, 
a = ^ ^ ^ (22) 



tr(R:,i^R: 



- 1 

-la!R_i ■ 



p = ^ — - " ' (23) 



7 



tr ( (R_, + M,)-^ + N,)) - trfRl,!^- 



1 2 



(24) 



Similar to (fTTl) . the numerator of « in (l22l) is not positive. If the bargaining set S is not empty, 
then Tj > . Thus, a is not positive. Similarily, /3 in (|23l) can be ignored as compared to 
a when the interference-plus-noise covariance matrices approach I. One can also see that 7 
in (l24l) is not positive. Therefore, when the interference-plus-noise covariance matrices R_j 
approach I, / "(t) < 0. 

The domain of the objective function in (fTSi) is {Qj|Qj ^ 0, tr(Qj) — < 0, % = 
1, . . . , M}, which is obviously convex. Therefore, the objective function in (fTSl) is concave. 

As a next step, let us prove the convexity of the constraint set. The constraint set of (fTSi) 
is identical to the bargain set ([8]). Specifically, it can be rewritten as 

^ = {Q,|Q, ^ 0, tr(Q,)-p,<0, ^ = l,...,M}n{Qi| -/i(Q) + jf^<0, ^ = 1,...,M} 
= 5in^2. (25) 
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It is easy to establish the convexity of the subset Si in (|25]) . but the convexity of the subset 
S2 is not obvious. Let us define h(t) = — /j(Q(t)) + I^^. Adopting the same methodology 
as the proof for the concavity of the utility function /^(Q), we have 

h"{t) = tr (^(R_, + M,)-^ + N,) (R_, + M,)-^ + N,) 

-Rll^R-^^l . (26) 



dt dt 

A similar result can be obtained that when the interference-plus-noise covariance matrices 
R_i approach I, the second term inside the trace operator in (l26l) can be ignored when 
compared to the first one. Thus, similar to (flTl) . if the interference-plus-noise covariance 
matrices R_j approach I, (l26l) can be rewritten as 

h"it) ^ tr((R_, + M,)-'(^ + N,^(R_, + M,r'(^^ + N, 

= tr ((Cf B,C,) (Cf B.Q)'') > (27) 

which implies that when the interference-plus-noise covariance matrices R_j approach I, 
h{t) is convex [23], i.e., the subset 5*2 in (1251) is convex. Consequently, the constraint set S 
in (1251) is convex as well. 

As we can see, the condition that the interference-plus-noise covariance matrices R_j 
approach I is sufficient for both the concavity of the objective function in (fTSi) and the 
convexity of its constraint set. Therefore, the Proposition 2 is proved. □ 

From Proposition 2, the following proposition can be obtained. 

Proposition 3: When the NB solution exists, the condition that the interference-plus-noise 
covariance matrices R^j approach I, is the sufficient condition for the uniqueness of the NB 
solution of MIMO interference systems. 

Proof: The interference-plus-noise covariance matrices R_j approaching I can ensure 
the convexity of the rate region (see Proposition 1) on which the NB is defined. Interestingly, 
it also guarantees the uniqueness of the NE in MIMO interference systems (see Proposition 1 
in [26]), i.e., it guarantees the IWF convergence. Moreover, Proposition 2 states that if the 
interference-plus-noise covariance matrices R_j approach I, then ([TSl) is a convex optimiza- 
tion problem, i.e., there exists at most one solution maximizing the objective function of 
(fTSl) [23]. Therefore, the requirement of the rate region convexity is the convexity of the 
optimization problem (fTSl) . Moreover, the existence of the solution for problem (fTSl) is the 
sufficient condition for the uniqueness of the NB solution. Applying Propositions 1 and 2 in 
this paper and Proposition 1 in [26], we complete the proof. □ 
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For proportional fairness, where the Nash product corresponds to the rate product, it is easy 
to infer from (fTSl) that the interference-plus-noise covariance matrices H i approaching I is 
also the sufficient condition for the concavity of ln(ln^ 2 Yliti -^«(Q))' which is equivalent to 
the concavity of the rate product logarithm for MIMO interference systems. 

2 ) Existence of NB Solution and Feasibility of NB Set: 

As stated in Section II, in MIMO interference systems, the NE is a steady and optimal state 
for selfish and competitive users competing with each other. However, the NE is not Pareto 
efficient, i.e., it is usually below the rate boundary. If all the users simultaneously deviate 
from the NE and agree to adopt different strategies (i.e., transmit covariance matrices), this 
may result in simultaneous rate improvement for all users. In this case, the NB solution will 
exist. 

The existence of the NB solution highly depends on the underlying MIMO interference 
channel, including its channel matrices, SNRs and INRs. In general, it is difficult to find a 
closed-form expression or quantitive description specifying when the NB solution exists. In 
the next section, we approach this problem by numerical studies and examine the impact of 
SNR and INR on the existence of the NB solution. 

Another important characteristic of the NB is its feasible bargaining set. Some general 
characteristics of the feasible utility set have been studied under an axiomatic framework in 
[28]. In what follows, we will analyze the structure of the transmit covariance matrices which 
form the feasible bargaining set of the NB in MIMO interference systems without employing 
time-sharing or orthogonal signaling. The transmit covariance matrix Qj (i = 1, - ■ ■ , M) of 
each user is Hermitian PSD. Thus, it can be decomposed as 

Q, = V,A,Vf, t = l,---,M (28) 

where Vj (i = 1, ■ ■ ■ , M) are unitary matrices of eigenvectors of Q, (i = 1, ■ ■ ■ , M) and 
Aj (z = 1, ■ ■ ■ , M) are the diagonal matrices with the corresponding eigenvalues. The transmit 
covariance matrices in the feasible set ([8]) can be considered as evolving from transmit 
covariance matrices of the NE (|7]). The evolution can be performed in three different ways: 
(i) changing the diagonal matrix (yUjI — Dj^^)+ of ([7]); (ii) modifying the unitary matrix U, 
of ©; (iii) changing both the diagonal and unitary matrices of (|7]). The feasible set S in (|25|) 
can be comprised of transmit covariance matrices of the following possible types: 

• Type I: The transmit covariance matrices Qj have the same unitary matrices Vj, but 
different diagonal matrices Aj from Q*. This is also referred to as stream control [17], i.e.. 
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the number of independent input data streams foreach user is limited. It can be interpreted as 
selectively shutting down some transmission antennas, which creates more interference 
to other users compared to the desired signal power generated for their target receivers. 

• Type II: The transmit covariance matrices have different unitary matrices Vj, but the 
same diagonal matrices Aj as Q*. This corresponds to the power control when all the 
effective users are the same, but their emitting directions are different from those of the NE. 

• Type III: Both the unitary matrices Vj and diagonal matrices of the transmit covari- 
ance matrices Qj for the NB are different from those of the NE. This can be viewed as 
a combination of Type I and II of the optimal transmit covariance matrices for the NB. 

C. Interference Cancellation in 2-User Systems 

Based on Remark 1, if the interference is high, interference cancellation techniques out- 
perform the orthogonal signaling techniques for guaranteeing the convexity of the rate region 
in the sense that applying interference cancellation leads to a larger rate region and, corre- 
spondingly, admits an NB solution with larger user rates than the orthogonal signaling based 
solution. Moreover, it ensures the uniqueness of NB solution. 

As a special case, we propose a multi-stage interference cancellation technique for 2-user 
systemso This technique will be used in the numerical studies section where the 2-user case 
is primarily considered. It is worth noting that our objective is to demonstrate the advantages 
of the interference cancellation over orthogonal signaling for the multicriteria optimization 
via NB rather than considering practically appealing interference cancellation designs. 

To design our interference cancellation technique, it is required that the receivers of both 
users cooperate with each other and perform multi-stage interference cancellation as follows. 
First, the receivers of users 1 and 2 decode their signals of interest and pass them to the 
receivers of users 2 and 1, respectively. Using ([T]), the decoded signals from receivers of 
users 1 and 2 in the first stage of interference cancellation can be written as: 



:HnVi 



- ^ -H^iy, (29) 



where x.^^ {i = 1,2) are the decoded signals at the receivers of users 1 and 2 in the 
first stage of the interference cancellation procedure. Second, after receiving the decoded 



''Note that the generalization to the multi-user case is straightforward, but it is omitted here for brevity. 
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signals from another user, the receivers of users 1 and 2 reconstruct the interference signals 
from the decoded signals, which can be expressed as ^/fj^^^^l2^i^^ and y^ry2iH2iXi^^ for 
users 1 and 2, respectively. Finally, each receiver subtracts the reconstructed interfering signal 
from its originally received signal so to obtain the following signals used in the second stage 
of the interference cancellation procedure 

y? = yi - Vrh2H.i2^2^ 

- f VpTHn - /^^HisH^-s^Hsi) Xi + m - . /^HisH^^s (30) 

y?^ = y2 - Vr^^2ixS^^ 

= (^Vp^H22 - y^^H2iHr/Hi2^ X2 + n2 - y^H2iHr>i (31) 

where y^-^^ (i = 1, 2) are the input signals at the second stage of interference cancellation. 
The procedure is repeated in the second stage to obtain the signals 

yf ^ = yi - v^Hi2 (^Vp^22 - y^^H2iHri^Hi2) yf 

= VPiHiiXi+ni-07i;Hi2^VP2H22 - ^J^^^li.2i^ii^v2j 

X ^n2 - y^H2iHri'ni^ (32) 

y? = y2 - V^H2i (^VpIHh - sj^^'i^i2^22^2^ yf ^ 

= -\/p2-^22X2 + n2 — -072lH21 ^y^Hn — .^^^^^Hi2H22^H21^ 

X ^ni-y^Hi2H2-2^n2^ (33) 

which are the input signals in the third stage of the procedure. The decoding, reconstruction, 
and subtraction are terminated after a certain stage or after a stopping criteria on the value of 
the remaining error is satisfied. Note that it is guaranteed that at each stage the new estimate 
of a signal is at least not worse than the estimate in the previous stage [29]. 

IV. Numerical Studies 

A. Convexity of the rate region 

Proposition 1 in Section EI provides the sufficient condition guaranteeing the convexity of 
the rate region for MIMO interference systems. Unlike Proposition 1 where the interference- 
plus-noise covariance matrices and then INRs are interpreted in a qualitative manner, in 
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this subsection we investigate the impact of the INRs on the convexity of the rate region 
quantitively via numerical studies. Consider a 2-user Rayleigh fading MIMO interference 
system with parameters a = 1 and Nt = Nr = 2. We examine the probability that the rate 
region is convex for different values of SNR/INR. It can be seen in Fig. 1 that when the SNR 
is 0.1, 5 and 50, the SNR/INR should be at least 10, 15 and 20dB, respectively, to ensure 
that the rate region is convex. Two trends can be observed in Fig. 1: (i) as the SNR/INR 
increases, i.e., as the INR becomes decreasingly small while the SNR is fixed, the probability 
that the rate region is convex increases; (ii) as the SNR increases, the SNR/INR needs to 
increase as well to retain the probability that the rate region is convex. 

When the Pareto rate region is nonconvex, we can resort to either mixed strategies, e.g., 
time-sharing or orthogonal signaling, or interference cancellation techniques to construct a 
convex rate region for the NB. Fig. 2 gives an example of these convexifying schemes and 
their resulting convex rate regions for a 2-user MIMO interference system with A^^ = N,,- = 2, 
Hii = diag(1.8, 1), H22 = diag(l, 1.8), H12 = H21 = diag(l, 1), = p2 = 20 , r]i2 = 7721 = 
14 and pi = P2 = 1- 

Fig. 2 shows the maximum achievable rates for both users in the cases when (i) pure strate- 
gies are used, i.e., neither time-sharing nor orthogonal signaling is adopteqj, (ii) time-sharing 
is employed, (iii) orthogonal signaling is employed, i.e., FDM and TDM are used, (iv) the 
multi-stage interference cancellation proposed in Subsection IIII-CI is applied. Specifically, a 
2-stage interference cancellation scheme given in (l32l) and (l33l) is considered. Its resulting 
rate boundary and the new NB solution is named as the IC boundary and NB solution 
after IC, respectively, in the figure. The NB solution which can be obtained graphically 
from the Nash curvqj and the NE for this MIMO system are also shown in Fig. 2. It 
can be seen from the figure that the orthogonal signaling schemes, especially the TDM, 
result in significant rate loss. The time-sharing signaling has a larger rate region than TDM 
and FDM. While, the multi-stage interference cancellation leads to the largest rate region. 
This example shows the potentials and advantages of interference cancellation techniques 
over the orthogonal signaling and time-sharing schemes in terms of guaranteeing a larger 

^For the sake of brevity, we name the resuUing rate region/boundary as the pure-strategy rate region/boundary in the 
sequel. 

*If the rates for users 1 and 2 at NE are denoted as and li"^, the Nash curve can be algebraically expressed 

as (y — l2^){x — Ii^) — max(Nash product). The NB solution can be then interpreted as the intersection of the rate 
boundary and the Nash curve as shown in Fig. 2. 
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convex rate region and, correspondingly, an NB solution with larger user rates. Moreover, 
the interference cancellation-based approach is more desirable than the time-sharing signaling 
due to the fact that the corresponding rates are actually achievable, while the rates promised 
by the time-sharing-based approach are achievable only in average. Note that practically 
instantaneous rates are more important than average rates and so the desired instantaneous 
rate must be achievable. 

B. Fairness of the NB 

NB is a bargaining approach which balances the individual fairness and system-level 
efficiency. The fairness of the NB is guaranteed by the fact that its resultant utility for 
each user is not less than that of the NE, which is considered as a relatively fair approach 
due to its selfish and competitive nature. In this subsection, we study the fairness of the NB 
solution and compare it with other bargaining approaches. The following bargaining solutions 
from [30] are taken as the benchmark for comparison: 

• Egalitarian solution: It is an absolutely fair solution with identical rate for each user. 

• Kalai-Smorodinsky (K-S) solution: It results in utilities proportional to their maximal 
achievable rates. 

• Utilitarian solution: This solution maximizes the sum rate of all users. 

• Proportional solution: It maximizes the rate product of all users. 

We compare the fairness of these bargaining solutions with the NB solution in terms of Jain's 
fairness index (JFI) [31] having the Egalitarian solution as the optimal solution. Recently, 
the JFI has been also used to compare difference scheduling algorithms in MIMO broadcast 
channels [32]. JFI J for an n-user system is obtained by Jain's equation [31] 

J = ^^fe^ (34) 

where Xi = Ti/T*, {Ti, Ts, • • • , T„} and {T*, T*, ■■■ , T*} are the measured and optimal 
utility vectors for the n-user system, respectively. It rates the fairness of the system and 
ranges from 1/n (worst case) to 1 (best case). 

To compare the fairness of the NB solution with that of the above mentioned bargaining 
solutions, we consider a 2-user MIMO interference system as an example with Nt — N^ — 2, 
Hn = diag(1.8, 1), H22 = diag(0.5, 0.6), H12 = diag(l, 1), H21 = diag(0.3, 0.3), pi = P2 = 
20 , r]i2 — 7721 = 14 and pi — P2 — Fig. 3 shows the pure-strategy rate boundary, time- 
sharing rate boundary and different bargaining solutions for this MIMO interference system. 
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These bargaining solutions are demonstrated graphically and geometrically in Fig. 3. For 
example, the NB solution can be interpreted as the intersection of the rate boundary and 
the Nash curve. Similarily, the Egalitarian solution, K-S solution. Utilitarian solution and 
Proportional solution are the intersections of the rate boundary and the following curves 
y = X, y/x = max(/2)/max(/i), y + x = max(/i + I2), yx = max(/i ■ I2), respectively. 
JFIs obtained from (|34|) for the NE, NB solution, K-S solution. Utilitarian solution and 
Proportional solution are 0.9925, 0.9685, 0.8853, 0.8960 and 0.8960, respectively. Note that 
in this example, the Utilitarian solution is identical to the Proportional solution. It can be seen 
that the JFI of the NB solution is only slightly smaller than that of the NE but much larger 
than all the other solutions. Compared with the NE, the NB improves the sum rate of this 
MIMO system by 18% at the price of only compromising its JFI by 0.9925-0.9685=0.024. 
These results confirm that the NB is an effective approach well balancing the indivisual 
fairness and system performance of interference systems. 

C. Feasibility of NB Set 

To exemplify the feasibility of NB set discussed in Section III-B, we first consider a 2- 
user MIMO interference system with the same setup as Fig. 2. In this example, the transmit 
covariance matrices of the NE obtained by the IWF are Q^^ = diag(0.25, 0.75) and Q^^ = 
diag(0.75, 0.25). The transmit covariance matrices resulting in the NB are Qf^ = diag(l, 0) 
and Q^^ = diag(0, 1). As we can see, the transmit covariance matrices Qf^ (i = 1,2) are 
of Type I as mentioned in the previous section, with the same unitary matrices but different 
diagonal matrices from Qf^. 

As a second example, consider a 2-user MIMO interference system with 

^2 2\ /l 0' 

Nt = Nr = 2, Hn = H22 = , H12 = H21 

V V v \ 

pi = p2 = 20, 7]i2 = f]2i = 60 and pi = P2 = 1- Fig. 4 demonstrates the pure-strategy rate 
boundary, the NE and NB solutions. It can also be seen from Fig. 4 that the rate region is 
convex, so there is no need to employ signaling. The transmit covariance matrices of the NE 
in this example are 

' 0.50 0.50 \ / -0.71 0.71 \ / \ / -0.71 0.71 
0.50 0.50 / \ 0.71 0.71 / \ 1 / \ 0.71 0.71 
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Whereas, the transmit covariance matrices leading to the NB solution are 

0.02 0.15\ /-0.99 O.ieA (oo\ f -0.99 0.16 
0.15 0.98 / I 0.16 0.99 / \ 1 / I 0.16 0.99 

In this scenario, the transmit covariance matrices Qf^ have the structure of Type II, i.e., they 
have the same diagonal matrices but different unitary matrices. Alternatively, the transmit 
covariance matrix can be interpreted as preceding. The corresponding preceding matrices for 
its NB and NE are 

piVB^ /O 0-16\ /O 0.71 

I 0.99 y ' ' \ 0.71 
As the last example, consider another 2-user MIMO interference system with 

1 o' 
^0 1^ 

P1—P2 — 20, 7712 = 7721 = 16 and pi—p2 — l. Fig. 5 shows the time-sharing rate boundary, 
pure-strategy rate boundary, the NE and NB solutions. It can be seen from Fig. 5 that the 
rate region is nonconvex in this MIMO system, but its NB solution still lies on the pure- 
strategy rate boundary. The transmit covariance matrices of the NE in this MIMO interference 
system are / 



/2 1^ 




(10) 




Uii 


, H22 = 




, H12 — H21 






i02i 





0.80 0.40 \ / -0.89 -0.45 \ M \ / -0.89 -0.45 
0.40 0.20 / \ -0.45 0.89 / \ / \ -0.45 0.89 



Q 



. 0.22 -0.14 \ / -0.23 0.97 
-0.14 0.78 / V 0.97 0.23 




V 



-0.23 0.97 
0.97 0.23 



and the transmit covariance matrices leading to the NB solution are 

0.90 0.29\ 1-0.95 -0.3l\ 1 1 o\ ( -0.95 -0.31 
0.29 0.10 / I -0.31 0.95 / \ / \ -0.31 0.95 



Qi 



NB 



0.21 -0.40 \ / -0.45 0.89 \ 10\ / -0.45 0.89 



\NB _ 

-0.40 0.79 / \ 0.89 0.45 / \ / \ 0.89 0.45 




- 

As we can see, Qf® and Q^^ have the structure of Types II and III, respectively, i.e., both 
the unitary and diagonal matrices for Q^^ are different from those of the corresponding NE. 
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D. The SNR and INR Impact on the Existence of the NB Solution 

Figs. 6 and 7 depict the SNR and INR impact on the NE and NB solutions, respectively. 
The channel realization is the same as that used for obtaining Fig. 4. It can be seen from 
Fig. 6 that the NB solution provides better rates as the SNR increases. Another phenomenon 
is that at low and high SNRs the NE coincides with the NB solution, i.e., the NE is the 
optimal solution. This is due to the fact that the NE lies on the rate boundary at these SNRs. 
In this example, the NE is optimal when SNR is less than 15 or larger than 230. In Fig. 7, 
the NB solution remains unchanged in the INR range where the NB solution exists. This 
is due to the fact that the NB in this case adopts stream control, i.e., only one stream is 
transmitting for each user, which eventually converts the MIMO interference system into an 
interference-free system. 

Fig. 8 presents the existence of the NB solution at different values of the SNR and INR. 
The symbol ' x ' represents the fact that the NB solution exists at the corresponding SNR and 
INR, i.e., the rate of the users in NE can be further improved by simultaneously changing their 
transmit covariance matrices from the NE. As mentioned in Section III-B, for systems with 
different channel matrices, the SNR and INR impact on the existence of their NB solutions 
may vary. In this MIMO interference system, the NB solution exists only when the INR is 
smaller than the SNR. From Figs. 7 and 8, it can also be seen that the NB solution does not 
exist when the INR is too small. The reasons behind this phenomenon are as follows: 

• When the INR is very small compared to the SNR, if stream control is adopted, the 
negative contribution to the rate from the loss in the desired signal component will be 
larger than the positive contribution from the decrease in interference. It is preferable 
for both users to compete selfishly with each other rather than adopting stream control. 
Therefore, in this very small INR region which is depicted as the selfish competition 
region in Fig. 8, the NE and NB solutions are identical. 

• As the INR increases and exceeds a certain value, the negative contribution to the rate 
from the loss in desired signal component becomes smaller than the positive contribution 
from the decrease in interference after adopting stream control, but the interference is 
still not strong enough to lead the IWF to converge to the stream control solution. In 
this case, the NB ends up with the stream control solution but the NE does not, and the 
NB solution exists. 

• As the INR continues increasing, the interference becomes stronger, which stimulates 
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the IWF to approach the stream control solution. When the INR reaches or exceeds a 
certain larger value, the NE converges to the stream control solution. In this region with 
large INR (we name it as the stream control region in Fig. 8), the NE coincides with 
the NB solution again. 
Another phenomenon can be observed from Figs. 6 and 8 is that the NB solution does not 

exist when the SNR is relatively small or large. This phenomenon can be explained similarly 

as the previous one: 

• When the SNR is small compared to the INR, the IWF is driven to the stream control 
solution which is also the optimal solution in that region, i.e., the NE and NB solutions 
are the same, both of them choose stream control. 

• As the SNR increases and exceeds a certain threshold, the IWF is driven to deviate from 
stream control and approach selfish competition, but the optimal (NB) solution in this 
area is still stream control. In this region, the NB solution exists and outperforms the 
NE solution which adopts selfish competition. 

• As the SNR continues increasing, the NB solution begins to approach selfish competition 
since the interference becomes less significant compared to the desired signal. When the 
SNR reaches a certain larger value, the NB solution converges to selfish competition 
which is same as the NE solution. 

E. Uniqueness of the NB Solution 

As a natural extension to our numerical studies of the existence of the NB solution, we also 
examine the uniqueness of the NB solution by calculating the value of /" (t) in (|2T]) when the 
NB solution exists. It can be seen from (EB-dal) that / (t) is a function of the underlying 
channel (channel matrices, SNRs and INRs). Meanwhile, it is also a function of t, (Zi, Z2) 
and (Xi, X2) in the convex combination (flOl) . Let the channel matrices be the same as those 
used to obtain Fig. 2. We set t = 0.5, and choose transmit covariance matrices (Zi,Z2) 
and (Xi, X2) such that their corresponding rates would be in the maximum achievable rate 
boundary, their rates for user 1 are I^^ and I^^ll, respectively. With this setup. Fig. 9 depicts 
the value of /"(t)|t=o.5. It can be seen that in this instance the value of f" {t) is negative in 
the region where the NB exists. We can infer from Section III-B that when the NB solution 
exists and the value of f"{t) is not positive for any t, (Zi,Z2) and (Xi,X2), then its NB 
solution is unique. 
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V. Special Cases 



In the previous section, we were concerned with characterizing the Pareto rate region and 
the NB over MIMO interference systems. Here, we further emphasize that some findings and 
results regarding the convexity of the rate region, the uniqueness, existence and the feasibility 
of the NB set have much broader validity than just in the MIMO case. Particularly, these 
results also hold for MIMO and SISO interference systems. In fact, the NB over MIMO 
interference systems unifies and generalizes that of MISO and SISO scenarios. In what 
follows, the applicability of the Pareto rate region and some NB characteristics is analyzed 
in the case of MISO and SISO. 

Remark 2: The sufficient conditions which guarantee the convexity of the rate region and 
the uniqueness of the NB solution for MIMO interference systems derived in Section III still 
hold in both MISO and SISO interference systems. 

The expression for the mutual information ([21) is a general one and is applicable for MIMO, 
MISO and SISO systems. The channel matrix in ^ converts from an A^^ x Nt matrix 
for MIMO to a 1 X A^t vector hjj and a complex value hij for MISO and SISO, respectively. 
The transmit covariance matrix Qj is an Nt x Nt matrix for both MIMO and MISO, and a 
non-negative real value Qi for SISO. More specifically, equation ^ can be written as 



The optimization problems Q and (fTSl) are still applicable to MISO and SISO systems. We 
further express the transmitted signal vector in ([T]) as x/ = FjSi, where Fj is an Nt x A^^ 
matrix which is known as precoding matrix in MIMO and Sj is the A^^ x 1 information symbol 
vector for user i. Without loss of generality, we assume that E[siS^] = I. Then, the transmit 
covariance matrix can be generalized as 



where Wi is the A'i x 1 beamforming vector for MISO transmitter i and Wi is the complex 
valued power control weight for SISO transmitter i. Thus, the optimization problem ^ 
can be interpreted as follows: the NB aims at designing the optimal precoding matrices, 
beamforming vectors and power control weights for all cooperative users under the power 




for MISO 



for SISO 



(35) 



E[FiFf] for MIMO 
Qi = < E[wiwf] for MISO 
E[\wi\^] for SISO 



(36) 
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constraints, with the objective to maximize the Nash product of the MIMO, MISO and SISO 
interference system, respectively. Therefore, sufficient conditions ensuring the convexity of 
the rate region and the uniqueness of the NB solution in MIMO interference systems derived 
in Section III-A and B can be extended to both MISO and SISO scenarios, i.e., Propostion 
1~3 are still valid for MISO and SISO interference systems. 

Remark 3: The analytical structure of the transmit covariance matrices leading to the 
NB in Section III-B and the numerical studies on existence of the NB solution in MIMO 
interference systems are still valid in MISO case. 

This is due to the fact that these two characteristics are analyzed based on the transmit 
covariance matrix Qj, which is an x N-^ Hermitian PSD matrix for both MIMO and 
MISO systems. However, as for a SISO interference system, the transmit covariance matrix 
Qj degenerates to a non-negative real value Qu so these results are not applicable to SISO 
interference systems. 

VI. Conclusions 

In this paper, the rate control problem in multi-user MIMO interference systems has been 
formulated as an MCO problem. The convexity of the Pareto rate region of this MCO problem 
has been studied. A sufficient condition which guarantees the convexity of the rate region 
has been derived. It is argued that the interference cancellation techniques are preferable 
for convexifying the Pareto rate region. Then, the MCO problem has been transformed into 
a single-objective optimization problem using NB. A variety of characteristics for the NB 
solution in MIMO interference systems such as the uniqueness and the existence of the 
NB solution and the feasibility of the NB set have been investigated. A sufficient condition 
ensuring the uniqueness of the NB solution in MIMO interference systems has also been 
derived. Moreover, a multi-stage interference cancellation scheme has been proposed to 
convexify the rate region of the corresponding MCO problem. It is shown that it leads to 
an NB solution with larger user rates. The convexity of the rate region, the effectiveness of 
the proposed interference cancellation scheme, the fairness of the NB solution, the impact 
of the SNR and INR on the NE and NB solutions, the existence of the NB solution, and 
its uniqueness (if exists) have also been demonstrated via numerical studies. Finally, the 
applicability of several NB characteristics to the MISO and SISO interference systems has 
been shown. 
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Fig. 3. Various bargaining solutions for a MIMO interference system. 
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Fig. 4. Rate region, NE and NB solution with Type II feasible NB set. 
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Fig. 5. Rate region, NE and NB solution with Type III feasible NB set. 
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Fig. 6. Impact of the SNR on the NE and NB solution (1NR=14). 
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Fig. 7. Impact of the INR on the NE and NB solution (SNR=20). 
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Fig. 8. Existence of the NB solution at different SNRs and INRs. 



